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AN EXTENSION OE DELSARTE’S METHOD. 
THE KISSING PROBLEM IN THREE AND 
EOUR DIMENSIONS 

Oleg R. Musin * 


1 Introduction 

The kissing number k(n) is the highest number of equal nonoverlapping spheres 
in R" that can touch another sphere of the same size. In three dimensions the 
kissing number problem is asking how many white billiard balls can kiss (touch) 
a black ball. 

The most symmetrical configuration, 12 billiard balls around another, is if 
the 12 balls are placed at positions corresponding to the vertices of a regular 
icosahedron concentric with the central ball. However, these 12 outer balls do 
not kiss each other and may all moved freely. So perhaps if you moved all of 
them to one side a 13th ball would possibly fit in? 

This problem was the subject of a famous discussion between Isaac Newton 
and David Gregory in 1694. (May 4, 1694; see interesting article m for details 
of this discussion.) It is commonly said that Newton believed the answer was 12 
balls, while Gregory thought that 13 might be possible. However, Bill Casselman 
1^ found some puzzling features in this story. 

This problem is often called the thirteen spheres problem. R. Hoppe m 
thought he had solved the problem in 1874. But, Thomas Hales CHI in 1994 
published analysis of Hoppe’s mistake (see also m)- Finally this problem was 
solved by Schiitte and van der Waerden in 1953 m A subsequent two-pages 
sketch of an elegant proof was given by Leech m in 1956. No much doubts 
that Leech’s proof is correct, but there are gaps in his exposition, many involved 
sophisticated spherical trigonometry. (Leech’s proof was presented in the first 
edition of the well known book by Aigner & Ziegler the authors removed 
this chapter from the second edition because a complete proof to include so 
much spherical trigonometry.) The thirteen spheres problem continues to be of 
interest, new proofs have been published in the last few years by Wu-Yi Hsiang 
|21| . Karoly Boroczky |2], and Kurt Anstreicher |2]. 
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Note that fc(4) > 24. Indeed, the unit sphere in centered at (0, 0,0, 0) 
has 24 unit spheres around it, centered at the points (±-\/2, ±-\/2,0, 0), with any 
choice of signs and any ordering of the coordinates. The convex hull of these 24 
points yields a famous 4-dimensional regular polytope - the “24-ceH”. Its facets 
are 24 regular octahedra. 

Coxeter proposed upper bounds on k{n) in 1963 for n = 4,5,6, 7, and 
8 these bounds were 26, 48, 85, 146, and 244, respectively. Coxeter’s bounds 
are based on the conjecture that equal size spherical caps on a sphere can 
be packed no denser than fc -|- 1 spherical caps on that simultaneously touch 
one another. Boroczky proved this conjecture in 1978 [S]. 

The main progress in the kissing number problem in high dimensions was 
in the end of 1970’s. Vladimir Levenshtein Ea, and independently Andrew 
Odlyzko and Neil Sloane [57|, Chap. 13] using Delsarte’s method in 1979 
proved that k{8) = 240, and fc(24) = 196560. This proof is surprisingly short, 
clean, and technically easier than all proofs in three dimensions. 

However, n = 8, 24 are the only dimensions in which this method gives a 
precise result. For other dimensions (for instance, n = 3,4) the upper bounds 
exceed the lower. In EZI the Delsarte method was applied in dimensions up to 
24 (see Table 1.5]). For comparison with the values of Coxeter’s bounds 
on k{n) for n = 4, 5,6, 7, and 8 this method gives 25, 46, 82, 140, and 240, 
respectively. (For n = 3 Coxeter’s and Delsarte’s methods only gave fc(3) ^ 13 
EHEZI ) Kabatiansky and Levenshtein have found an asymptotic upper bound 
20 . 401 »( 1 -Ho( 1 )) fQj. ^ bound 20-2075n(l+o(l)) 

in m 

Improvements in the upper bounds on kissing numbers (for n < 24) were 
rather weak during next years (cni Preface to Third Edition] gives a brief review 
and references). Arestov and Babenko |3] proved that the bound fc(4) ^ 25 
cannot be improved using Delsarte’s method. Hsiang m claims a proof of 
fc(4) = 24. His work has not received yet a positive peer review. 

If M unit spheres kiss the unit sphere in R", then the set of kissing points 
is an arrangement on the central sphere such that the (Euclidean) distance 
between any two points is at least 1. So the kissing number problem can be 
stated in other way: How many points can be placed on the surface of S"“^ so 
that the angular separation between any two points is at least 60° ? 

This leads to an important generalization: a finite subset X of S"“^ is called 
a spherical z-code if for every pair (x, y) of X the scalar product x-y < z. Spheri¬ 
cal codes have many applications. The main application outside mathematics is 
in the design of signals for data transmission and storage. There are interesting 
applications to the numerical evaluation of n-dimensional integrals [HU Chap.3]. 

The Delsarte method (also known in coding theory as Delsarte’s linear 
programming method, Delsarte’s scheme, polynomial method) is described in 
cni EH Let fit) be a real polynomial such that fit) < 0 for t G [—1,2;], the 
coefficients c^’s in the expansion of fit) in terms of Gegenbauer polynomials 
are nonnegative, and cq = 1. Then the maximal number of points in a 
spherical 2 :-code in is bounded by /(I). Suitable coefficients c/c’s can be 
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found by the linear programming method [Tni Chapters 9,13]. 

We found an extension of the Delsarte method in 2003 m (see details in 
|2ti|l. that allowed to prove the bound fc(4) < 25, i.e. fc(4) = 24. This extension 
yields also a proof fc(3) < 13. 

The first version of these proofs was relatively short, but used a numerical 
solution of some nonconvex optimization problems. Later on |2ti| these calcula¬ 
tions have been reduced to calculations of roots of polynomials in one variable. 
(This is not a big problem now, all computer algebra systems such as Maple, 
Mathematica, and Matlab can find roots. Also these calculations can be inde¬ 
pendently verified. If you have approximate values all roots of a polynomial, 
then you can check the existence of these roots by simple computations.) 

We present in this paper a new proof of the Newton-Gregory problem, an 
extension of Delsarte’s method, and a proof that fc(4) = 24. 


2 The thirteen spheres problem: a new proof 

Let us recall the definition of Legendre polynomials Pk (t) by recurrence formula: 

„ 3nl „ 2A;-1„ k-1 ^ 

Pq = ^, Pi = t, P2 = -t Pk = — t Pk-i - - — Pk-2; 

2 2 K k 

1 

or equivalently Pfc(t) = , ,, ~ 1)^ (Rodrigues’ formula). 

2^ k\ df^ 

Lemma 1. Let X = {xi, X 2 , ■ ■ ■, Xn} be any finite subset of the unit sphere 
in R^. By fiij = dist(a:i,Xj) we denote the spherical (angular) distance between 
Xi and Xj. Then 

n n 

0 . 

i=i j=i 

Let 

, 2431 n 1287 7 18333 , 343 4 83 o 213 n t 1 
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Lemma 2. Suppose X = {xi,X2, ■■■ ,Xn} C . Then 


i=i j=i 




Lemma 3. Suppose X = {xi, X 2 ,..., x„} is a subset o/S^ such that the angular 
separation fiij between any two distinct points xt^xj is at least 60°. Then 


= X! X! /(cos(())ij)) < 13n. 

i=i i=i 
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Theorem 1. k{3) = 12. 

Proof. Suppose X is a kissing arrangement on with n = fc(3). Then X is 
satisfying the assumptions in Lemmas 2, 3. Therefore, nf ^ S{X) < 13n. From 
this follows n < 13, i.e. n ^ 12. From other side we have fc(3) ^ 12, then 
n = fc(3) = 12. □ 

We need the only one fact from spherical trigonometry, namely the law of 
cosines: 

cos (f) = cos 9i cos 02 + sin 6 i sin 62 cos tp, 

where for spherical triangle ABC the angular lengths of its sides are 0i, 6 * 2 , </> and 
the angle between AB,AC is ip (Fig. 1). li ip = 90°, then cost/) = cos^i cos 02 
(spherical Pythagorean theorem). 


A 




-1 -0.5 0 0.5 0.8 

Fig. 2. The graph of the function /(t) 


Proof of Lemma 1. 

This lemma easily follows from Schoenberg’s theorem m for Gegenbauer 
polynomials. Note that Pk = G\ . For completeness we give a proof of Lemma 
1 here. In this proof we are using original Schoenberg’s proof that based on the 
addition theorem for Gegenbauer polynomials.^ 

The addition theorem for Legendre polynomials was discovered by Laplace 
and Legendre in 1782-1785: 


k 

Pfe(cos 01 cos 02 + sin 01 sin 02 cos(/j) = ^ Cm,/c P™(cos 0 i)P™(cos 02 ) cosrmp 

m—0 

= Pfc(cos0i)Pfc(cos02) + 2 ^ P™(cos0i)P™(cos02) cosrmp, 

^' (k + m)\ 

m—1 ^ ' 

^Pfender and Ziegler 1 ^ give a proof as a simple consequence of the addition theorem for 
spherical harmonics. This theorem is not so elementary. The addition theorem for Legendre 
polynomials can be proven by elementary algebraic calculations. 


4 








where 


(See details in |S1E].) 

Proof. Let X = C and Xi has spherical (polar) coordinates 

(0i, ipi). Then from the law of cosines we have: 

cos4 >i^j = cos 9i cos 9j + sin 9i sin 9j cos^, ipij = i.pi — tpj, 

which yields 

k 

Pk{cos(j)ij) = X! y] <^m,kPJP'{cos9i)PJP'{cos9j) cosm(pij 

i,j i,j m—0 

— ^ ^ Cm,k ^ ^ COS Ul^i^j , Um.i ~ ^k (cOS 9i^. 

m i,j 

Let US prove that for any real ui^... ^Un 

UiUj cos rrupij ^ 0. 

i,3 

Pick n vectors yi,... ,yn in with coordinates yi = (cos rrupi, sin rrvpi). If 
y = uiyi + ... + M„y„, then 

< y,y >= I |y| 1^ = y] UiUj cos rrvpij ^ 0. 

ij 

This inequality and the inequalities Cm,k > 0 complete our proof. □ 


Proof of Lemma 2. 

Proof. The expansion of / in terms of Pk is 
9 

/ = CkPk = Po + 1 . 6 Pi + 3.48P2 + 1.65P3 + 1.96^4 + O.IP 5 + O. 32 P 9 . 

k=0 

We have cq = 1, ^ 0, fc = 1,2,..., 9. Using Lemma 1 we get 

9 n n n n 

s{x) = Y,ckY.Y. Pfe(cos(0i,j)) > EE cqPo = 
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Proof of Lemma 3. 


Proof. 1. The polynomial f{t) satisfies the following properties (see Fig.2): 

(i) f{t) is a monotone decreasing function on the interval [—1, —to]; 

(ii) fit) <0 for t e (-< 0 ,1/2]; 
where /(—to) = 0, to « 0.5907. 

These properties hold because /(t) has the only one root —to on [—1,1/2], 
and there are no zeros of the derivative f'{f) (8th degree polynomial) on [—1, —to]. 


Let SiiX) := ^ /(cos)/)^^)), then SiX) = ^ SiiX). From this follows 

i=l 

if S'i(X) < 13 for i = l,2,...,n, then S{X) < ISn. 

We obviously have (fi^i = 0, so /(cos(/i,i) = /(I). Note that our assumption 
on X {(fij ^ 60°, i ^ j) yields coscfij < 1/2. Therefore, coscfij lies in the 
interval [-1,1/2]. Since (ii), if coscpij € [—to, 1/2], then f{cos(j)ij) < 0. Let 
'^(*) ■= {/ ■ coscfij G [—1,—to)}. We obtain 

S^iX) < T,iX) := /(I) + ^ /(cos/).,,). (1) 

Let 9o = arccosto,6lo « 53.794°. Then j G J{i) iff /).,,■ > 180° — 9o, i.e. 
9j < 9q, where 9j = 180° — /).,,-. In other words all Xij, j G J(t) lie inside the 
circle of center eo and radius 9o, where eo = —Xi is the antipodal point to x.. 

2. Let us consider on points eo, j/i,..., j/m such that 


= dist{yi,yj) > 60° for all i j, dist(eo,yi) < 6»o for 1 ^ t < m. (2) 

Denote by y the highest value of m such that the constraints in (2) define a 
non-empty set of points yi,... ,ym- 

Suppose 0 ^ m ^ /i and Y = {j/i,..., y^} satisfies (2). Let 

H{Y) = H{yi, ...,y.^):= /(!) + /(_ cos6»i)-f.. . + fi-cos9m), 9i = dist(eo,2/i) 

hjn ■= maxi7(y), Umax := max{/io, hi,.h^}. 

It is clear that Ti{X) ^ hm, where m = |J(i)|. From (1) it follows that 
SiiX) ^ hm- Thus, if we prove that hmax < 13, then we prove Lemma 3. 

3. Now we prove that /i ^ 4. 

Suppose Y = {yi ,..., ym} C satisfies (2). If eo is the North pole and yi has 
polar coordinates (0., ipt), then from the law of cosines we have: 


cos (fij = cos 9i cos 9j + sin 9i sin 9j cos(v3i — ^Pj). 
From (2) we have cos/).,,- ^ 1/2, then 


cos((/?i 


Pj) < 


1/2 — cos 9i cos 9j 
sin 9i sin 9j 


( 3 ) 
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Note that arccos((l/2 — to)/(l ~ ^o)) ~ 76.582° > 72°. Then m ^ 4 because 
no more than four points can lie in an unit circle with the minimum angular 
separation between any two points greater than 72°. 

4. Now we have to prove that hmax = max{/io, hi, / 12 , / 13 , hi} < 13. 

We obviously have ho = /(I) = 10.11 < 13. 

From (i) follows that /(— cos0) is a monotone decreasing function in 9 on 
[ 0 , 6 * 0 ]. Then for m = 1 : H{yi) = /(I) + /(—cos0i) attains its maximum at 
01 = 0 , 

^ 1 =/(I)+ /(-!) = 12.88 <13 

5. Let us consider for m = 2, 3,4 an optimal arrangement {eo, 1 / 1 ,..., ym} 
in that gives maximum of HiY) = hm- Note that for optimal arrangement 
points yk cannot be shifted towards eo because in this case H{Y) increases. 

For m = 2 this yields: eo € yij/ 2 , and dist(j/i, 1 / 2 ) = 60°. If eo ^ yiy 2 , then 
whole arc 7 / 11/2 can be shifted to ep. Also if dist(j/i, 7 / 2 ) > 60°, then 7/1 (and 7 / 2 ) 
can be shifted to cq. 

For TO = 3 we prove that A 3 = 7 /i 7 / 2 j /3 is a spherical regular triangle with 
edge length 60°. As above, eo G A 3 , otherwise whole triangle can be shifted 
to eo. Suppose dist( 7 /i, t/^) > 60°, i = 2, 3, then dist( 7 /i,eo) can be decreased. 
From this follows that for any yi at least one of the distances {dist( 7 /i, 7 /^)} is 
equal to 60°. Therefore, at least two sides of A 3 (say 7 / 17/2 and 7 / 17 / 3 ) have length 
60°. Also dist( 7 / 2 , 2 / 3 ) = 60°, conversely 7/3 (or 7 / 2 , if ep G 7 / 17 / 3 ) can be rotated 
about 7/1 by a small angle towards eo (Fig.3). 

When TO = 4 first we prove that A 4 = 7 / 17 / 27 / 37/4 is a convex quadrangle. 
Conversely, we may assume that 7/4 G 7 / 17 / 22 / 3 . 

The great circle that is orthogonal to the arc 602/4 divides into two hemi¬ 
spheres: Hi and H 2 . Suppose eo G Hi, then at least one 2/i (say 2 / 3 ) belongs H 2 
(Fig.4). So the angle Zeo 2 / 42/3 greater than 90°, then (again from the law of 
cosines) dist( 7 / 3 , ep) > dist( 7 / 3 , 2 / 4 ). Thus, 

03 = dist( 7 / 3 , 60 ) > dist( 2 / 3 , 2 / 4 ) ^ 60° > 0 o — a contradiction. 

Arguing as for to = 3 it is easy to prove that A 4 is a spherical equilateral 
quadrangle (rhomb) with edge length 60°. 
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6 . Now we introduce the function where f/' G [60°, 29o]: 

Fi(V’) := max Fi(6», V-) =/(-cos6») +/(-cos(V’- 0)). 

'IP/2^0^0q 

So if dist{yi,yj) = then 

f {-cos 0i) + f{-cos Oj) ^ Fi{ip). (4) 


Therefore, 


H{yi, y2) ^h2 = /(I) + i^i(60°) « 12.8749 < 13. 

7. When m = 4, A 4 is a spherical rhomb. Let di = dist(j/i, j/ 3 ), and 
d 2 = dist(j/ 2 , 2 / 4 ), then cos((ii/2) cos((i2/2) = 1/2 (Pythagorean theorem, the 
diagonals of A 4 are orthogonal). So if p{s) := 2 arccos[l/(2 cos(s/2))], then 
di = p{d2), d2 = p{di), p{90°) = 90°. 

Suppose di ^ d 2 - Since 0i ^ Oq, d 2 ^ 20q, then 
p(2^q) ^ di ^ 90° ^ ^2 ^ 20q. 

Now we consider two cases: 

1) p{20o) < di < 77°, and 2) 77° ^ di < 90°. 

1) Fl{^p) is a monotone decreasing function in ip. Then (4) implies 
/(-cosdi) + f{-cos 03 ) < Fi{p{20o)), /(-cos 6 * 2 ) +/(-cosd 4 ) < Fi(p(77°)), 
so then 


H{Y) < /(I) + Fiip{20o)) + Fiip{77°)) « 12.9171 < 13. 
2) In this case we have 

H{Y) ^ /(I) + i^i(77°) + Fi(90°) « 12.9182 < 13. 


Thus, /14 < 13. 

8 . Our last step is to show that /13 < 13.^ 

^For given ip, the value F\{ip) can be find as the maximum of the 9th degree polynomial 
11 ( 5 ) = _Fi( 0 , ' 0 ), s = cos (0 — '?/?/2), on the interval [cos( 0 o ~ !]■ 

^More detailed analysis shows ^3 12.8721, /i 4 12.4849. 




Since A 3 is a regular triangle, H{Y) = /(I) + /(—cos0i) + /(—cos 02 ) + 
/(— cos 6 * 3 ) is a symmetric function in 0 ^, so we can consider only the case 

^1 ^ ^ ^3 ^ ^ 0 - 

In this case Rq ^ O 3 ^ Oq, where Rq = arccos -\/2/3 « 35.2644°. (Note that 
the circumradius of A 3 equals Rq.) 

Let yc is the center of A 3 . Denote by u the angle Zcoysyc- Then (see Fig.5) 
cos 01 = cos 60° cos 03 + sin 60° sin 63 cos (i?o — u), 

cos 02 = cos 60° cos 03 + sin 60° sin O 3 cos {Rq + u), 

where Zyiysy^ = Zy 2 y 3 Vc = Ro, 0 < m < mq = arccos(cot 63/ VS) - Rq 
(if u = uq , then 02 = 03). 

For fixed 63 = 'll;, H{yi,y2) becomes the polynomial of degree 9 in s = cosm . 
Denote by F2{ip) the maximum of this polynomial on the interval [cosuo, 1]. 
Let 

{V-i,..., V'e} = {Ro, 38°, 41°, 44°, 48°, 0o}. 

It’s clear that F 2 (' 0 ) is a monotone increasing function in ip on [i?o,0o]- From 
other side, f{—cos ip) is a monotone decreasing function in ip. Therefore for 
03 e [ipi,ipi+i] we have 

H{Y) = H{yi,y2) + /(- cos03) < Wi := F2{ipi+i) + /(- cosipi). 


Since, 


{wi,..., ws} ~ {12.9425,12.9648,12.9508,12.9606,12.9519}, 

we get h 3 < maxjryi} < 13. 

Thus, hm < 13 for all m as required. 


□ 


3 Delsarte’s method 

Let X = {a;i, a; 2 ,..., xm} be any finite subset of the unit sphere C R", 
S”“^ = {x : X £ R", x-x =\\x\\^ = 1}. From here on we will speak of a: G 
alternatively of points in 8"“^ or of vectors in R”. 

By (pij we denote the spherical (angular) distance between Xi, Xj. It is clear 
that for any real numbers ui,U 2 , ■ ■ ■, um the relation 

11 UiXi I P = cos (pijUiUj > 0 

id 

holds, or equivalently the Gram matrix T{X) is positive semidefinite, where 

T^X) — (tij), tij — COS (pij — Xi ' Xj. 

Schoenberg m extended this property to Gegenbauer (ultraspherical) poly¬ 
nomials of Uj. He proved that if gij = G^^\tij), then the matrix (g^) is 
positive semidefinite. Schoenberg proved also that the converse holds: if f{t) is 
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a real polynomial and for any finite X C the matrix is positive 

semidefinite, then / is a sum of with nonnegative coefficients. 

Let us recall the definition of Gegenbauer polynomials. Suppose (t) be 
the polynomials defined by the expansion 

OO 


Then the polynomials = G^"^(t)/G^"^(l) are called Gegenbauer or ul- 

traspherical polynomials. (So the normalization of G^”^ is determined by the 
condition G^"^(l) = 1 .) 

Also the Gegenbauer polynomials g["^ can be defined by recurrence formula: 


g(") = 1, 


G^”^ = t, 


Gt^ = 


(2fc + n-4)tGt\-(fc-l)Gt^, 
fc + n — 3 


They are orthogonal on the interval [—1,1] with respect to the weight func¬ 
tion p(t) = (1 — t 2 )(Ti- 3)/2 ^ggg details in [Rl UniIT?)ll^ 1 . In the case n = 3, G^"^ 

are Legendre polynomials Pfe, and G^^^ are Ghebyshev polynomials of the second 
kind (but with a different normalization than usual, C/fc(l) = 1 ), 


Gf{t) = Uu(t) 


sin((A: -f 1 ) 0 ) 
(fc -I- 1) sin0 ’ 


t = cos 0 , 


fc = 0,1,2,... 


For instance, Uq = 1, Ui = t, U 2 = (4t^ — l)/3, U 3 = 2t^ — t, 

U 4 = - 12^ + l)/5, ..., Gg = (256t9 - h\2t^ -f 336^® - 5t)/5. 

Let us now prove the bound of Delsarte’s method. If a matrix {gij) is positive 
semidefinite, then for any real ut the inequality '^gijUiUj ^ 0 holds, and then 

for Mi = 1, we have ^ 0- Therefore, for gij = G^^\tij), we obtain 

i,j 

M M 

(3-1) 

i=i j=i 

Suppose 

f{t) = coG^^\t) + ... + CdG^^\t), where Co ^ 0,..., Cd ^ 0. (3.2) 


Let S{X) = E E fitij)- Using (3.1), we get 
^ j 


d M M MM 

^ EEco^o"^(^b) = COM2. (3.3) 

fc—0 i—lj—l i—lj—1 

Let X = {xq, ... ,xm} C S”“^ be a spherical z-code, i.e. for all i ^ j, 
tij = cos(pij = Xi ■ Xj ^ z, i.e. tij G [—l,z] (but tu = 1). Suppose f{t) ^ 0 for 
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te[-l,z], then S'(X) = M/(l) + 2/(<i2) + ... + 2/(tM-iM) < M/(l). If we 
combine this with (3.2), then for cq > 0 we get 

M ^ (3.4) 

Co 

The inequality (3.4) play a crucial role in the Delsarte method (see details in 
HEiisiiiniiisiisiiiaEii^). If 2 = 1/2 and Co = 1, then (3.4) implies k{n) ^ 
/(I). In 121 EH Levenshtein, Odlyzko and Sloane have found the polynomials 
f {t) such that /(I) = 240, when n = 8; and /(I) = 196 560, when n = 24. 
Then k{S) ^ 240, fc(24) ^ 196 560. When n = 8, 24, there exist sphere packings 
{Eg, and Leech lattices) with these kissing numbers. Thus A:(8) = 240 and 
fc(24) = 196 560. When u = 4, a polynomial / of degree 9 with /(I) = 25.5585... 
was found in (22]. This implies 24 ^ A:(4) ^ 25. 


4 An extension of Delsarte’s method. 

Let us now generalize the Delsarte bound M ^ /(l)/co. 

Definition. Let f{t) be any function on the interval [—1,1]. Consider on 
points j/o, yi, • ■ ■, 2/m such that 

2 /j -Vj ^ z for all i ^ j, f{yo ■ yi) > 0 for 1 ^ < m. (4.1) 

Denote by p, = y,{n, z, /) the highest value of m such that the constraints in 
(4.1) define a non-empty set of points (yoj • • • j 2/m)- 
Suppose 0 ^ m ^ yt. Let 

H{Y) = H{yo;yi,...,ym) ■= f{^) + fivo ■ Vi) + ■ ■ ■ + fivo -ym), 

hm ■■= max{iL(T)}, hmax ■= max {ho, hi ,... 

Remark, hmax depends on n, z, and /. Throughout this paper it is clear what 
/, n, and z are; so we denote by hmax the value hmax{n, z, /). 

Theorem 2. Suppose X C 8"“^ is a spherical z-code, |X| = M, and 
f{t) = coG^\t) -b ... -I- CdG^^\t), where cq > 0, ci > 0,..., > 0. Then 

M ^ IhUfzl. = i_ uiax{ho, hi,..., h^}. 

Co Co 

Proof. Since / satisfies (3.2), then (3.3) yields 

Six) ^ cqM^. 

Let J{i) := {j : f{xi • xj) > 0, j ^ z}, and X{i) = {xj : j € Then 

M 

Si{X) ^'^f{xi ■ Xj) ^ f{l) + ^ f{xi- Xj) = H{xi;X{i)) ^ hmax, 

1=1 j&Jii) 
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so then 


M 

sm = i: S'*(X) ^ Mhmax- 
i=l 

We have cqM^ ^ S{X) ^ Mhmax, i-e. cqM < hmax as required. □ 


Note that /iq = /(!)■ If f{t) ^ 0 for all t G [—1, z], then for a z-code X we 
have fJ, = 0, i.e. hmax = ho = /(!)• Therefore, this theorem yields the Delsarte 
bound M ^ /(l)/co. 

The problem of evaluating of hmax in general case looks even more compli¬ 
cated than the upper bound problem for spherical z-codes. It is not clear how 
to find /i? Here we consider this problem only for a very restrictive class of 
functions f{t): f{t) ^0 for t G [—to,z\, to > ^ > 0 . 

Let us denote by A(k, to) the maximal number of points in a spherical s-code 
n C of minimal angle uj, cosco = s. (Note that H(n,60°) is the kissing 
number k{n).) 


Theorem 3. Suppose Y = {yi, . ■ ■ ,ym} is a spherical z-code in S"“^, and 
points Pi lie inside the sphere of center eo and radius 9o, where to = cos 6*0 ^ z. 
Then 


m ^ n — I, 


arccos 


z-tl \ 

l-tl)- 


Proof. We have cj)i,j = dist(yi, yj) ^ S = arccos z for i j; 

6 i = arccos(eo • Pi) ^ 9o for 1 ^ z ^ to; and 6 *o ^ 6 . 

Let n be the projection of Y onto equator from pole eo- Denote by 'jij 
the distances between points of H in S"“^. Then from the law of cosines and 
the inequality cos (pij Gi z, we get 


cos 7i J- 


cos (j)ij — cos Oi cos 0 j 
sin 6 i sin 6 j 




z — cos 9i cos 9j 
sin 9i sin 9j 


Let Q{a) 


z — cos a cos f3 
sin a sin (3 ’ 


then Q\a) 


cos P — z cos a 
sin^ a sin P 


From this follows, if 0 < a,/3 < then cos/3 > z (because 9o ^ 6); so then 
Q'{a) p 0, and Qia) p Q{9o). Therefore, 


cos7ij 




z — cos 9i cos t 
sin 9i sin 0 , 


< 


z — cos (to 
sin^ 9o 


z — t 
1 - t 


2 

0 

2 

0 


that complete our proof. □ 

Corollary 1. Suppose f{f) < 0 for t G [—to,z], toP z p 0, then 


y{n,z, 



n — 1 , arccos 


z-to \ 

l-tp- 
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Proof. The assumption on / yields f{yo ■ Ui) > 0 only if 
0 i = dist(eo,yi) < 0q = arccosto, 

where cq = —yo is the antipodal point to yo. Therefore, this set of points 
{Co, 2 / 1 ,, 2/m} satisfies the assumptions in Theorem 3. □ 

The next claim will be applied to prove that A: (4) = 24. 

Corollary 2. Suppose f{t) ^ 0 for t G [—to, 1/2], to ^ 0.6058, then 
M = m(4, 1/2,/) ^ 6. 

Proof. Note that for tg > 0.6058, arccos[(l/2 — to)/(l — tp)] > 77.87°. So 
Corollary 1 implies /t(4,1/2, /) ^ A(3, 77.87°). 

Denote by (pk{M) the largest angular separation that can be attained in 
a spherical code on containing M points. In three dimensions the best 

codes and the values (p 3 {M) presently known for M ^ 12 and M = 24 (see 
Ha El EOl) For instance, Schiitte and van der Waerden m proved that 
1 / 53 ( 5 ) = 1 / 53 ( 6 ) = 90° and (^ 3 ( 7 ) ~ 77.86954° (cos/53(7) = cot40° cot 80°). 

Since 77.87° X/ 53 (7), then A(3,77.87°) < 7, i.e. /x < 6. □ 

Corollary 1 shows that if to is close enough to 1, then /x is small enough. 
Then one gets relatively small - dimensional optimization problems for compu¬ 
tation of numbers hm for small n. If additionally f{f) is a monotone decreasing 
function on [—1,—to], then these problems can be reduced to low-dimensional 
optimization problems of a type that can be treated numerically. 


5 Optimal sets for monotonic functions 

In this section we consider f{t) that satisfies the monotonicity assumption: 
f{t) is a monotone decreasing function on the interval [—1, —to], 

fft) ^0 for t G [-to,z], to> z^O (*) 

Consider on 8"“^ points 2 / 0 ,2/i, ■ ■ ■ ,ym that satisfy (4.1). Denote by 0k for 
fc > 0 the distance between x/fc and Cq, where cq = —yo is the antipodal point to 
2 / 0 . Then 2/0 • 2/fc = — cos 0k, and H(Y) is represented in the form: 

HiY) = /(I) -b /(- cos6*i) -b ... -b /(- cos6i^). (5.1) 

A subset C of S"“^ is called (spherical) convex if it contains, with every two 
nonantipodal points, the small arc of the great circle containing them. If, in 
addition, C does not contain antipodal points, then C is called strongly convex. 
The closure of a convex set is convex and is the intersection of closed hemispheres 
(see details in USD- If a subset Z of S" ^ lies in a hemisphere, then the convex 
hull of Z is well defined, and is the intersection of all convex sets containing Z. 
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Suppose f{t) satisfies (*), then Qm = {yi, ■ ■ ■, Vm} lies in the hemisphere of 
center cq. Denote by the convex hull of Qm in = convQm- 

Now we consider an optimal arrangement of Qm for H. Let 6 = arccosz, 
(jjij = dist(j/i, yj), N{Qm) = number of (j)ij =6 {yi ■ y^ = z). 

Definition We say that Qm is optimal if H{Y) = hm- If optimal Qm is not 
unique up to isometry, then we call Qm as optimal if it has maximal N{Qm)- 

The function f{t) is monotone decreasing on [—1,—tg]- By (5.1) it follows 
that the function H{Y) increases whenever 9k decreases. This means that for 
an optimal Qm no yk € Qm can be shifted towards eg. 

That yields 

eg G Am (5-2) 

because in the converse case whole Qm can be shifted to eg. 

From this follows that for m = 1, eg = j/i. Thus 

hi = f{l) + f{-l). 

It was proved in Section 2 that for m = 2 : dist(?/i, 2 / 2 ) = <5, thus 
/12 =/(I) + max {(/(— COS0) +/(—cos(d — 0))}, 0g = arccostg. 

6/2^0^0q 

It was also proved that A 3 is a spherical regular triangle with edge length 6. 
Using similar arguments it’s not hard to prove that for n > 3, A 4 is a spherical 
regular tetrahedra with edge length 6. 

Let Am, m ^ n, is a spherical regular simplex with edge length 6, and 

^771 = {y - y & Am, yyk^to, 1 ^ fc ^ m}. 

Note that flm is a convex set in Let 

Hm{y) = /(I) + f{-y ■yi) + ... + f{-y ■ ym). 

Then hm is the maximum of Hm{y) on ^Im- 

hm = max Hm{y), ^m C Am C 2 ^ TO ^ min(n,/i). (5-3) 

veAm 

When n > m any yk G Qm is a vertex of Am- In other words, no yk that lies 
inside Am- In fact, that has been proved in Section 2 (see 5, Fig. 4). 

In the first version of the paper m has been claimed that 
for optimal Qm with m > n, for any yk € Qm there are at least n — 1 distinct 
points in Qm o-t the distance of S from yk- 

However, Eiichi Bannai and Makoto Tagami found some gaps in our exposition. 
Most of them are related to “degenerated” configurations. In this paper we need 
only the case n = 4, to = 5. For this case they verified each step of our proof, 
considered all “degenerated” configurations, and finally gave clean and detailed 
proof. I wish to thank Eiichi Bannai and Makoto Tagami for this work. Now 
this claim in general case can be considered only as conjecture. 

^For m ^ n, is a spherical regular simplex with edge length 5. In this paper we need 

just cases m = 3,4. 
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6 An algorithm for computation suitable poly¬ 
nomials f{t) 

In this section is presented an algorithm for computation “optimal” ® polynomi¬ 
als / such that f{t) is a monotone decreasing function on the interval [—1, —to]j 
and f{t) ^0 for t G [—to, 2 ], to > z ^ 0. This algorithm based on our 
knowledge about optimal arrangement of points yi for given m. Coefficients Ck 
can be found via discretization and linear programming; such method had been 
employed already by Odlyzko and Sloane m for the same purpose. 

d . 

Let us have a polynomial / represented in the form f{t) = l+J2 {t). 

k=l 

We have the following constraints for /: (Cl) Cfe ^ 0, 1 ^ k ^ d; 

(C2) f{a) > f{b) for —1 ^ a < 6 < —to] (C3) /(t) < 0 for —to ^t ^ z. 

When m ^ n, hm = maxiLm(y), y € ^m- We do not know y where Hm 
attains its maximum, so for evaluation of hm let us use yc — the center of A^- 
All vertices yk of are at the distance of Rm from y^ where 
cos Rm = ^/{1 + {m— l)z)/m. 

When m = 2n — 2, Am presumably is a regular (n — l)-dimensional cross¬ 
polytope. (It is not proven yet.) In this case cos Rm = V^. 

Let/„ = {1, ..., n} lj{2n — 2}, m G In, bm = —cos Rm, whence 
Hm{yc) = /(I) + mf {bm)- If Fo is such that Hm{y) ^ E = Fq + /(I), then 
(C4) f{bm) ^ Fo/m, m & In- A polynomial / that satisfies (C1-C4) and gives 
the minimal E (note that E = Fo + 1 + ci + ... + Cd = Fo + f {!) will become a 
lower estimate of hmax) can be found by the following 

Algorithm. 

Input: n, z, to, d, N. 

Output: Cl, ...,Cd, Fo, E. 

First replace (C2) and (C3) by a finite set of inequalities at the points 
aj=-l + ej, Os^js^N, e = {l + z)/N: 

Second use linear programming to find Fq, ci,..., Cd so as to minimize 

d 

E — 1 = Fo + ^k subject to the constraints 

fc=i 


Cfe ^ 0, 1 ^ fc d; 








uj e [-1,-to]; 


d d 

l + ^CfeG^”^(aj) ^ 0 , ajG[-to,z]] 1+ '^CkG^,^\bm) ^ Fo/m, mein- 

k^l k^l 

Let us note again that E = max Hm{yc) ^ h^ax here, and that E = hmax 

m£ln 

only if hmax = HmoiVc) for some mo G In- 

^Open problem: is it true that for given to, d this algorithm defines / with minimal hmax^ 
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7 On calculations of hm for m ^ n 

Here we explain how to solve the optimization problem (5.3). Let C 8 ™“^ is 
a spherical regular simplex with edge length 6 = arccos z; yi, i = 1 ,..., to, are 
the vertices of A^; U = y ■ Vi = cosOi ^ to = cos^o; to > z; f{t) is a monotone 
decreasing function on the interval [— 1 ,—tg]; hm is the maximum of Hm{y) 
subject to the constraints U ^ to; Hm{y) = /(I) + f{-y ■yi) + ... + f{-y ■ ym)- 

The first method. 

Hm{y) is a symmetric function in the variables 0i,..., Om- Then we can consider 
this problem only on the domain K = {y ■. 6 m ^ ^ 62 ^ hi]. Note that A is 

a spherical simplex. Let us consider a barycentric triangulation of this simplex 
such that the diameter of any simplex ai of this triangulation is not exceed e. 

It is easy to prove that for any yk, y ■ yk attains its maximum on ai at some 
vertex of Ui. Denote this vertex by yk,i. Let I = {i ■. yi^i • yi ^ to}- So for i € / 
we have 

f{-yk,i ■ yk) = max{/(-j/ • yk)}, 

yetTi 

then 

m 

hm ^ • yk)Y 

* k^l 

That yields a very simple method for calculation of hm- For / from Section 
9 this method gives /13 « 24.8345, ^4 « 24.818. 

The second method. 

For m ^ n the values hm can be calculated another way. We are using here that 
/(t) = /o + /it + . • . + fdt‘^ is a polynomial. The first method is technically easier 
then the second one. However, the second method doesn’t assume that / is a 
monotone decreasing function on [— 1 ,—tg], and it can be applied to functions 
without monotonicity assumption. 

Let us consider Hmiy) as the symmetric polynomial Fm(ti,... ,tm) in the 
variables U = y ■ yi : Fm{ti, ...,tm) = /(I) + /(-ti) + ... + /(-tm)- Denote by 
Sfc = Sk{ti,..., tm) the power sum + ... + Then 

Fmih, ...,tm)= , Srf) = /(I) + to/o - /l Si + . . . + (-Ij'^/dSrf. 

From the fact that A^ is a spherical regular simplex follows 

® P-1) 

Any symmetric polynomial in to variables can be expressed as a polynomial 
of si,..., Sm- Therefore, in the case k > m the power sum Sk is Rk{si ,..., Sm). 
Combining this with (7.1), we get 

^m(si, Cr(Si), S 3 , ... , Sd) = $m(si, S 3 

1 • • • 1 ^m')- 

Therefore, we have 

hm — max $77^(51, S3,..., Sttt.) , (si, S3, . . . , Stt^,) € hJm C R. , 
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where Dm is the domain in dehned by the constraints U ^ to and (7.1). 

Let us show now how to determine Dm for m > 2. The equation (7.1) defines 
the ellipsoid E : S 2 = cr(si) in space {ti,..., tm}- Then si = ti +.. . + tm attains 
its maximum on E at the point with ti = t 2 = ■ ■ ■ = tm, and si achieves its 
minimum on Ef^{ti ^ to} at the point with t 2 = ■■■ = tm = to- From this 
follows ^ Si ^ W 2 , where 

ViP-to)iP- f l + {m-2)z 

wi = - |-(m-l)to, p= -^-, 

p m — 1 

W2 = a/to (to — 1) 2; + TO. 

The equation si = u) gives the hyperplane, and the equation S 2 = cr(w) 
gives the (to — l)-sphere in space: {(ti,..., tm)}. Denote by S{uj) the (to — 2)- 
sphere that is the intersection of these hyperplane and sphere. Let lk{oj) be the 
minimum of Sk on S{u!) ^ ^o}) and Vk{oj) is its maximum. Now we have 

hm = max max... max<I)m(si, S 3 ,..., Sm), where 

Si S3 Sm 


101^31^102, Ik(si) ^ Sk < Vk(si), fc = 3, ...,TO. 

For the polynomial / from Section 9 (and Section 2) we can give more details 
about calculations of hm for to = 3,4. 

Let us consider the case to = 3 with d = 9. In this case ^^(sa) = <I)3(w, S3) 
is a polynomial of degree 3 in the variable S3. 

Lemma 4. Let f be a 9th degree polynomial f(t) = /o + fit + ... + /gt® such 
that fg >0, fe = fs = 0, cmd fr > — 15 / 9 / 7 . If Ff{s) <0 at s = hitv), then the 
function FLj(s) achieves its maximum on the interval [ 13 ( 0 ;), ^ 3 ( 0 ;)] at s = 

Proof. The expansion of sg in terms of s^s^s*, t + 2j + 3fc = 9, is 


S9 = gS3 


+ 'S3(3'5i + ^ZSl) + S3(-S2 - + R{si,S2)- 


The coefficient of s§si in S 7 equals 7/9. Thus 


Fb(s) = -s^ fg /9 - s^ (/g wct(w) + 2/9 w ^/3 - 7/7 w/ 9 ) + sR/w) + i?o(w). 


F^{s) is a cubic polynomial with negative coefficient of s^. Then F^(s) is a 
concave function for s > r, where r : F”{r) = 0. Therefore, if r < Ig^Lo), then 
F^{s) is a concave function on the interval [/ 3 (a;), U 3 (a;)]. r < Is^to) iff 

B^lo) := 3^3(0;) + + 9 a; cr(a;) > —Icofr/fg. 


This inequality holds for tg < —z < 0. Indeed, 

uj wi ^ 1 + 2z, a (oj) > I, ^ 3 ( 0 ;) > 0; 


so then 

B(co) > 15a; > -luofr/fg- 


The inequality Ff{ls{uj)) ^ 0 implies that F^{s) is a decreasing function on the 
interval [I 3 (a;),ua(a;)]. □ 
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The polynomial / from Section 9 satisfies the assumptions in this lemma. 
Then $ 3 (w, s) attains its maximum at the point s = h^Lo), i.e. at the point with 
ti = t 2 ^ ts, or with ti ^ t 2 > ts = tf). If = ^2 ^ ta, then p(w) = <I> 3 (w, 
is a polynomial in to. This polynomial is a decreasing function in the variable to 
on the interval fa ^ to- Therefore, p{lo) achieves its maximum on this interval 
at the point with ts = to. The calculations show that for / from Section 9 
ho = maxp(w) « 24.8345, when 63 = Oq, 9i = 02 ^ 30.0715°. 

Corollary 3. Let f be the polynomial from Section 9, then ho « 24.8345. 


Consider the function F^{so, S 4 ) = ^ 4 ( 10 , S 3 , S 4 ) on S{uj). Let qi G S{uj) and 
qi : ti = t 2 > to = t 4 , q 2 '■ h = t 2 = to > ^4, and qo : h > t 2 = to = t4- 

Lemma 5. Let f be a 9th degree polynomial f(t) = ^fith ///g > 0 and 
fe = fs = Oj then the function TIj(s 3 ,S 4 ) achieves its maximum on S(uj) with 
u) > 1 at one of the points (soiqi), S 4 {qi)), i = 1,2, 3. 

Proof. The expansion of sg in terms of s^^s^s^s^ is 


Sg — Yg'® 4 '®l + n '®3 ~ q + 7-545351 + -S4S2S;i — 753823^ — 77535 ^ +i?(si, S2). 
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The coefficient of s|si in S 7 equals 0. We have /e = /s = 0, then 
Fajiso,S 4 ) = -/gSg + ... = -/g(s|/9 - slui^/S) + ... Therefore, 


_ 9^ifo(s3,S4) 
d^so 


r ^2 23 , 

-79(353 - 



- 53). 


If Fc^{so,S 4 ) has its maximum on S{u)) at the point x, and x is not a critical 
point of S 3 on S{lo), then F 33 < 0. From other side, for all ti G [0,1] and 
Si = w > 1 we have S 3 ^ w < so then F 33 > 0. The function S 3 on S{uj) 
(up to permutation of labels) has critical points a,t qt, i = 1 , 2 ,3. □ 


Corollary 4. Let f be the polynomial from Section 9, then /14 « 24.818. 

Proof. By direct calculations it can be shown that 

TL(s 3 (gi), 54 (^ 1 )) > Fuj{so{qi), S4{qi)) for i = 2,3. Then Lemma 5 implies 
h 4 = maxp(w), where p{uj) = Fuj{so{qi), S 4 {qi)) = 4 ) 4 (w, so{qi), S 4 (gi)). 

The polynomial p{to) attains its maximum ft ,4 « 24.818 at the point with 
01 = 6»2 « 30.2310°, 00 = 04 - 51.6765°. □ 


8 On calculations of in four dimensions 

Let us consider the case n = 4, m = 5. For simplicity here we consider only the 
case z = 1/2. Then d = 60° and 0o = arccostg < 60°. 

Denote by Fg the graph of the edges of A 5 with length 60° , where Qs is an 
optimal set. The degree of any vertex of Fs is not less than 3 (see Section 5). 
This implies that at least one vertex of Fs has degree 4. Indeed, if all vertices 
of Fg are of degree 3, then the sum of the degrees equals 15, i.e. is not an even 
number. There exists only one type of Fs with these conditions (Fig. 6 ). 
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For fixed dist (?/25 Vi) = ct, Q 5 is uniquely defined up to isometry. Therefore, 
we have the 1-parametric family A 5 (a) on S^. If dist(j/ 3 , j/ 5 ) = {3, then 

2cosacos/3 +COSO + cos/3 = 0 (8.1) 

The equation (8.1) defines the function (3 = X{a). Then a = A(/3), A(90°) = 90°. 
For all i we have dist(j/i, cq) ^ Oq, then 

dist( 2 /i, yj) ^ dist(?/i, eo) + dist(yj, cq) ^ 20o- 

Suppose a ^ /3, then (8.1) and the inequality /3 < 29q yield 

ao ^ ck ^ 90° ^ /3 ^ 26*0, c(o ■= max{60°, A(26*o)}. 


Let 

H 5 iy, a) = /(I) + fi-y ■ 2/1 (a)) + • • ■ + f{-y ■ y^ia))- 

Then 

/15 = maxjiJo)?/, a)}, y G S^, y-yk{a)^to, 1 ^ fc < 5, oq ^ a ^ 90° (8.2) 

y,a 

We have four-dimensional optimization problem (8.2). Our first approach 
for this problem was to apply numerical methods m For the polynomial / 
from Section 9 this optimization problem was solved numerically by using the 
Nelder-Mead simplex method: H^{y,a) achieves its maximum 24.6856 at 
a = 60° and y with 6»i « 42.1569°, 6*2 = 6*4 « 32.3025°, 63 = 05 = 0 o- 
(The similar approach for the case n = 4, m = 6 gives the 3-parametric family 
A6(a,/3, 7 ), and for / from Section 9: « 22.5205.) 

Note that (8.2) is a nonconvex constrained optimization problem. In this 
case, the Nelder-Mead simplex method and other local improvements methods 
cannot guarantee finding a global optimum. It’s possible (using estimations 
of derivatives) to organize computational process in such way that it gives a 
global optimum. However, such kind solutions are very hard to verify and some 
mathematicians do not accept such kind proofs. Fortunately, an estimation of 
/i 5 can be reduced to discrete optimization problems. 
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Let dist(yi,y) = ip, and = f{-y ■ Vi) + f{-y • %)• It’s clear 

that for fixed attains its maximum at some point that lies in the 

great 2-sphere that contains yi,yi,yj. Now we introduce the function 
Suppose yiyiyj is a spherical triangle in with dist{iji,yi) = dist(j/i,) = 
60°, dist( 7 /i, yj) = 7 , denote by F{ip, 7 ) the maximum of ^i,j{y, V') on subject 
to the constraints y • yk ^ to, k = i,j- Then ^ij{y,'ip) ^ F('ip,j), so then 
d>2,4(j/,V’) < $ 3 . 5 ( 2 /,'*/') < F{'>p,f3). Thus 

^5 ^ /(I) + /(- cos V*) + F{p:, a) + Fpip, \{a)). (8.3) 

Let uq < ai < ... < ak < Ofe+i = 90°. It’s easy to see that F{\p,'^) is a 
monotone decreasing function in 7 . That implies for a G [ai,ai+i] : 

F{'ip,a) ^ F{'ip,ai), F{'ip,X{a)) ^ F^ip, X{ai+i)). Therefore, from (8.3) follows 

ho ^ /(I) + /(- cos-ip) + max {F{ip, ai) + F^-tp, A(ai+i))}. (8.4) 

O^i^k 

Note that (8.4) to reduce the dimension of the optimization problem (8.1) 
from 4 to 2. It is not too hard to solve this problem in general case. However, 
the polynomial / from Section 9 satisfies an additional assumptions that allowed 
to find a weak bound on /15 even more easier. 

Let us briefly explain how to check the following assumptions for /: 

1 ) $jj ( 7 /, ■;/)) achieves its maximum at one of the ends of the arc uiptp,^), where 

w(V’,7) := { 2 / : 2/ G dist(j/i,y) =■;/', yyeP^to, 1 = i,j}; 

2 ) Fpip,^) is a monotone increasing function in tp. 

For given 7 (7 = dist{yi,yj)) and ip the function ^ij{y,'ip) becomes a poly¬ 
nomial p(s) of degree d on [so,l], where s = cosu, u = Zj/ij/ij/c, and yc is 
the center of yiyiyj (see Section 2, 8 ). Then 1) holds iff p'{s) has no roots on 
(sq, 1 ), either if s : p'{s) = 0 , then p"{s) > 0 . 

Using 1) it’s easy to check 2). For the polynomial f{t) from Section 9 if 
7 > 62.41°, then p{s) achieves its maximum at s = sq (be. dist(j/j, y) = 9q), so 
it’s clear that 2) holds. From other side if 7 < 69.34°, then the arc uj{ip,j) lies 
inside the triangle yiytyj, therefore F{ip,^) increases whenever ip increases. 

Note that 1) gives us the explicit expression for F{ip,^) = max(p(so),p(l))- 
For fixed 7 and ip ^ipi from 2) follows F{ip,j) ^ F{ipi,^). 

Denote by ipL(i), '4’uii) tb® lower and upper bounds on ip that defined by 
the constraints a G a^+i], y • Pq ^ to, q = 1 ,... ,b. 

Let ipL(i) = ipi,o < ‘ipi,! < ■ ■ ■ < ipi,e < ipi,e+i = i’uii)- Recall that /(- cosip) is 
a monotone decreasing function in ip. Then 2) and (8.4) yield 

ho ^ /(I) + max max }, (8.5) 


where 

Ri,] = /(- cosipij) + F(ipij+i,ai)} + F{ipij+i,X{ai+i)). 

It’s very easy to apply this method. Here we need just to calculate the 
matrix (Ri.j) and the maximal value of its entries gives the bound on / 15 . For 
®F('i/,,60°) = F2('|/>) — /(I) (see Section 2, 8, Fig. 5). 
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/ from Section 9 and to ~ 0.60794, Oq = arccosto ~ 52.5588°, /(—to) = 0, this 
method gives the bound < 24.8434.^ 

Now we show how to find an upper bound on /ig. Let {cq, ?/i, ..., yo} G S^, 
H{yi, ...,ye) = /(I) + /(- cos6»i) + ... + /(- cosdg), where Oi = dist(eo,yi). 
Suppose 9i ^ 02 ^ ^ Qq. Now we prove that 9q ^ 45°. That can be 

proven as Corollary 2 (Section 4). Conversely, all 9i < 45°. In this case 
to* = COS06 > l/v^, and uj = arccos [(1/2 — t§,j)/(l — tg^)] > 90°. But if 
u > 90°, then A(3,w) ^ 4 (see ^01 El) - a contradiction. (In fact we proved 
that 05 > 45° also.) 

Let us consider two cases: (i) 0o ^ 06 > 50° (ii) 50° ^ 06 > 45°. 

(i) H{yi,...,yo) = iJ( 2 /i,..., j/ 5 ) +/(-cos06). We have 

H{yi,...,yo) ^ < 24.8434, /(-cos06) < /(-cos50°) « 0.0906, 

then iJ(j/i,... , 2 / 5 ) < 24.934. 

(ii) In this case all 9i ^ 50°. Therefore, we can apply (8.5) for 0o = 50°. This 
method gives 0 . 5 ( 50 °) < 23.9181, then H{yi,... ,yo) < 05(50°) < 23.9181, so 
then 


H{yi ,..., ?/6) < 23.9181 + /(- cos45°) « 23.9181 + 0.4533 = 24.3714 
Thus 

06 < max{24.934, 24.3714} = 24.934 


9 k{4) = 24 

For n = 4, 2 ; = cos 60° = 1/2 we apply this extension of Delsarte’s method with 

/(t) = 53.76t® - 107.52t^ + 70.56t® + 16.384t‘‘ - 9.8321^ - 4.128t^ - 0.434t - 0.016 

The expansion of / in terms of Uk = is 

f = Uo + 2 Ui + 6.121/2 + 3.4841/3 + 5.121/4 + 1.051/9 

The polynomial / has two roots on [—1,1]: ti = —to, tg « 0.60794, t 2 = 1/2, 
/(t) ^ 0 for t G [—to, 1/2], and / is a monotone decreasing function on the 
interval [— 1,—to]. The last property holds because there are no zeros of the 
derivative /'(t) on [— 1,—to]. Therefore, / satisfies (*) for 2 ; = 1/2. 

Remark. The polynomial / was found by using the algorithm in Section 6. 
This algorithm for n = 4, z = 1/2, d = 9, N = 2000, to = 0.6058 gives 
E « 24.7895. For the polynomial / the coefficients Ck were changed to “better 
looking” ones with E « 24.8644. 

achieves its maximum at a = 60°, 1 / ~ 30.9344°. Note that this bound exceeds the 
tight bound on ^ 24.6856 given by numerical methods. 
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We have to > 0.6058. Then Corollary 2 gives ^ < 6. Consider all m < 6. 
ho = /(I) = 18.774, h, = /(I) + /(-I) = 24.48. 
h 2 = /(I) + ^^max^^{(/(- cose) + /(- cos(60° - 0))} w 24.8644, 
where Oq = arccosto ~ 52.5588°. 

Note that h 2 can be calculated by the same method as in Section 2. Here 
^2 = /(I) + 2/(— cos30°) also. 

In Sections 7, 8 have been shown that 

ho « 24.8345, hi « 24.818, ho < 24.8434, ho < 24.934. 

Theorem 4. fc(4) = 24 

Proof. Let X be a spherical 1/2-code in with M = k{A) points. The poly¬ 
nomial / is such that hmax < 25, then combining this and Theorem 2, we get 
fc(4) ^ hmax < 25. Recall that k{A) ^ 24. Consequently, A:(4) = 24. □ 

10 Concluding remarks 

The algorithm in Section 6 can be applied to other dimensions and spherical 
z-codes. If to = 1, then the algorithm gives the Delsarte method. E is an 
estimation of hmax in this algorithm. 

Direct application of the method developed in this paper, presumably could 
lead to some improvements in the upper bounds on kissing numbers in dimen¬ 
sions 9, 10, 16, 17, 18 given in [23 Table 1.5]. (“Presumably” because the 
equality hmax = E is not proven yet.) 
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In 9 and 10 dimensions Table 1.5 gives: 

306 < k{9) ^ 380, 500 ^ A:(10) ^ 595. 

The algorithm gives: 

n = 9: deg f = 11, E = hi = 366.7822, to = 0.54; 
n = 10 : degf = 11, E = hi = 570.5240, to = 0.586. 

For these dimensions there is a good chance to prove that 
fc(9) < 366, A:(10) < 570. 

From the equality fc(3) = 12 follows :^3(13) < 60°. The method gives 
1 ^ 3 ( 13 ) < 59.4° (deg/ = 11). The lower bound on :^3(13) is 57.1367° [r7|. 
Therefore, we have 57.1367° < </53(13) < 59.4°. 

The method gives 1 ^ 4 ( 25 ) < 59.81°, Lpi{2A) < 60.5°. (This is theorem that 
can be proven by the same method as Theorem 4.) That improve the bounds: 

(/?4(25) < 60.79°, (/?4(24) < 61.65° El (cf. H); ipi{2A) < 61.47° H; 

(^ 4 ( 25 ) < 60.5°, (^ 4 ( 24 ) < 61.41° H. 

Now in these cases we have 

57.4988° < (^ 4 ( 25 ) < 59.81°, 60° ^ Lpi{2d) < 60.5°. 

For all cases that were considered {z < 0.6) this method gives better bounds 
than Fejes Toth’s bounds for m and Coxeter’s bounds for all (pn{M) 

El However, for n = 5,6,7 direct use of this generalization of the Delsarte 
method does not give better upper bounds on k{n) than the Delsarte method. 
It is an interesting problem to hnd better methods. 
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